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ABSTRACT 

This  investigation  is  a  study  of  the  problem  faced  by  a  circuit 
designer  who  in  the  analysis  and  synthesis  of  his  circuit  must  use 
parameters  which  vary  over  wide  ranges. 

A  technique  utilizing  a  digital  computer  resulting  in  design 
curves  is  developed.   A  second  technique  is  demonstrated  where  the 
same  design  curves  are  obtained  by  another  method  employing  an  analogue 
computer. 

Both  techniques  are  applied  to  a  neon  plot oconductor  circuit  and 
design  curves  are  obtained  with  each  method. 

The  writers  wish  to  express  their  appreciation  for  the  assistance 
and  encouragement  given  them  by  Professor  W.  Conley  Smith  of  the  U.  S. 
Naval  Postgraduate  School. 
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1.   Int roduet ion, 

Systew  reliability  is  the  keynote  in  this  dynamic  age  of  mechaniz- 
ation where  unfailing  performance  is  expected  and  required.  However, 
component  parts  in  a  systems  configuration  are  often  subject  to  manufactur- 
ing tolerances  which  may  vary  radically  from  the  norm,  yet  a  systems  de- 
signer has  to  deal  with  these  varying  components  and  produce  a  dependable 
system.  Dependability  can  be  measured  as  a  statistical  quantity  or  as  an 
absolute  requirement.  The  latter  definition  of  dependability  will  be  used 
in  this  paper  as  a  yardstick  in  measuring  the  success  of  a  particular  de- 
sign. The  problem  to  be  undertaken  in  this  thesis  is  the  synthesis  of  a 
reliable  system  consisting  of  widely  varying  components.  Reliability  as 
used  here  is  dependent  ©n  the  components  themselves  not  failing. 

Circuit  synthesis  is  generally  interpreted  as  the  formulation  of  sig- 
nal flow  charts  and  the  determination  of  component  values  for  the  realiza- 
tion of  a  given  transfer  or  imroittance  function.  This  paper  will  be  con- 
cerned with  only  the  latter  part  of  this  definition. 

The  specific  system  to  be  investigated  and  designed  consists  of  a 
neon  photoconductor  logic  circuit  which  is  a  basic  building  block  for  other 
logic  circuitry.  Justification  of  the  important  requirement  of  reliability 
in  this  circuit  is  indicated  by  its  use  as  a  basic  element  in  the  founda- 
tion of  an  entire  system.  Analysis  and  synthesis  of  the  circuit  will  be 
restricted  to  a  steady  state  investigation. 

Using  the  specific  neon  photoconductor  logic  circuit  as  an  example, 
a  general  technique  of  circuit  synthesis  will  be  proposed  in  which  the 
tools  of  circuit  analysis  and  inductive  reasoning  are  utilized  for  synthesiz- 
ing. Analysis  can  exist  without  synthesis  but  the  converse  of  this  state- 
ment is  not  true.  H.  S.  Scheffler  and  F.  R.  Terry  in  their  paper  entitled 
"Description  and  Comparison  of  Five  Computer  Methods  of  Circuit  Analysis" 

1 


describe  three  non-statiatieal  approaches  to  circuit  analysis.  Some  of 
their  concepts  in  circuit  analysis  are  adopted  and  extended  into  the  field 
of  circuit  synthesis.  For  example  the  criteria,  point  of  failure,  is  a 
basic  extreme  condition  applicable  to  a  circuit  problem  whether  the  investi- 
gation be  of  analysis  or  snythesis.  The  particular  logic  circuit  to  be 
investigated  will  be  synthesized  with  the  aid  of  both  a  digital  and  an 
electronic  analogue  computer.  Each  computer  has  advantages  over  the  other 
which  are  peculiar  to  the  methods  utilized  in  the  solution.  Successful 
circuit  synthesis  by  the  computers  produces  design  curves,  which  when 
properly  interpreted  can  be  an  invaluable  asset  to  the  circuit  designer. 


20  Description  of  Digital  Computer  Design  Technique,, 

Thif  ssetion  explains,  in  general  terms,  a  technique  useful  in  cir- 
cuit design  which  uses  to  full  advantage  the  high  speed  calculating  feature 
of  a  digital  e©mputere 

In  designing  a  circuit  to  meet  output  specifications  the  designer  is 
almost  invariably  faced  with  the  problem  of  using  components  In  the  cir- 
cuit which  do  not  have  constant  values,,  but  may  have  values  in  a  given 
range.  The  simplest  example  of  this  point  is  an  ordinary  fixed  resistance 
which  may  have  any  value  in  a  range  above  and  below  its  nominal  value* 
This  variation  is  due  to  manufacturing  differences.  A  rang*  of  variation 
may  also  be  due  to  environmental  changes  during  circuit  operation,  aging 
or  other  factors e 

Circuit  parameters  which  must  be  selected  from  within  a  known  range 
of  variation  will  be  referred  to  as  '"known"  parameters,,  Other  parameters, 
whose  value  can  be  chosen  by  the  designer  based  on  his  design  to  meet 
specifications  are  defined  as  "design"  parameters,, 

The  final  circuit  ©r  circuits  decided  upon  by  the  designer  must  meet 
the  specifications  for  any  random  set  of  the  "known"  parameters  -  each 
parameter  in  the  set  picked  arbitrarily  from  anywhere  in  its  range  ©f 
variation.  The  designer  must  choose  tha  design  parameters  for  his?  circuit 
so  that  when  "acceptable"  values  of  the  design  parameters  are  used  in  the 
circuit,  the  circuit  performance  will  meet  specifications  for  any  random 
set  of  the  known  parameters. 

In  simple  circuits  the  decision  as  to  acceptable  values  of  the  design 
paraneters  may  be  easy.  Consider  a  simple  voltig©  source  and  resistance 
circuit.  If  the  voltage  is  a  known  parameter  with  variation  from  5  to  10 
volts  and  th-2   specification  on  output  demands  a  current  of  fror.  3  t©  A 
amperes,  then  the  selection  of  acceptable  values  of  resistance^  the  design 


parameter 
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If  the  design  parameter,  R,  has  any  acceptable  value;  that  is  any 
value  from  1.25  to  3.33  ohtns,  the  circuit  will  meet  performance  specifica- 
tions for  any  value  of  the  known  parameter,  voltage,  in  its  range  of  5  to 
10  volts. 

The  example  shows  the  significance  of  the  terms  known  and  design 
parameters,  but  is  much  too  simple  to  demonstrate  the  technique  developed. 
In  the  general  case,  there  are  many  known  parameters,  each  with  a  range  of 
variation.   There  is  usually  a  number  of  output  specifications  to  be  met 
and  the  designer  is  at  liberty  to  pick  his  design  parameters  in  sets  of 
acceptable  values. 

Stated  most  simply,  the  technique  begins  by  first  writing  equations 
which  express  the  output  as  a  function  of  the  design  and  known  parameters. 
By  inspecting  these  equations  and  employing  boundary  conditions,  the  de- 
signer must  establish  a  range  of  investigation  for  each  of  the  design 
parameters.   This  range  is  similar  to  the  range  of  variation  of  the  known 
parameters,  but  differs  in  that  it  may  be  made  larger  or  smaller  during 
the  process  of  design  based  upon  the  designers  interpretation  of  the  in- 
itial answers. 

Once  a  range  has  been  established  for  each  parameter,  the  equation 
expressing  the  output  must  be  solved  for  each  value  in  the  range  of  each 
parameter.   The  values  of  design  parameters  which  yield  acceptable  output 
regardless  of  the  values  of  the  known  parameters  used  are  recorded  and 


plotted  us  acceptable  eowbinat ,!        :  design  parameters*  k  circuit 
built  with  an   acceptable  combination  of  these  desigr?  parameters  will  pro- 
vide output  which  meets  specifications  no  matter  what  values  of  known 
parameters  are  used* 

it 

Since  it  is  impossible  to  solve  the  circuit  equation  for  every  value 
of  every  design  parameter  and  every  known  parameter,  the  ranges  of  these 
parameters  must  be  divided  into  increments  and  the  equation  solved  for 
each  increment  of  each  parameter*,  The  selection  of  these  increments  in- 
troduces some  error*.  As  in  other  approximation  procedures,  smaller  in- 
crements introduce  less  error,  However,  smaller  increments  lead  to  more 
combinations  to  investigate  and  thus  many  more  calculations,   If  A  is 
the  number  of  increments  of  the  n  th  parameter  to  be  investigated,  the 
total  number  of  answers  to  be  computed  and  tested  is 

(A)  (a  ,..    //4j 

It  is  seem  that  any  reasonably  complicated  design  problem  with  a 
practical  number  of  increments  quickly  leads  to  an  extremely  large  number 
cf  calculations,.  This  problem  suggests  use  of  a  digital  computer  with 
its  rapid  calculating  ability  as  an  invaluable  toole 

In  some  cases  it  may  not  be  necessary  to  investigate  all  increments 
of  a  known  parameter.  This  exception  was  demonstrated  in  the  trivial  ex- 
ample stated  above  although  not  specifically  pointed  out.  The  range  of 
acceptable  resistance  was  determined  by  investigating  only  the  maximum 
and  minimum  values  of  the  known  parameter,  voltage.   Because  of  the  simple 
nature  of  the  circuit  it  is  obvious  that  a  value  of  the  design  parameter 
which  is  acceptable  at  both  the  maximum  and  minimum  values  of  the  known 
parameter  will  be  acceptable  for  any  value  between  and  only  the  end  points 
need  be  used.  The  characteristic  of  the  parameter  voltage  which  makes 
this  true  is  that  the  partial  derivative       is  always  of  the  same  sign. 


In   general  then,  if  the  partial  derivative  of  the  output  taken  with 
respect  to  a  known  parameter  is  always  of  the  same  sign,  for  all  values 
of  all  the  other  known  and  design  parameters,  only  the  maximum  and  mini- 
mum values  of  that  known  parameter  need  be  investigated,   since  if  the 
design  parameters  are  selected  to  be  acceptable  with  both  these  points, 
they  will  be  acceptable  for  any  value  of  this  parameter  in  its  range. 

Since  parameters  of  this  sort  mean  a  great  reduction  in  the  number 
of  calculations,  the  equations  for  output  should  be  partially  different- 
iated with  respect  to  each  known  parameter  in  turn  to  determine  if  deri\a- 
tives  of  constant  sign  exist.  The  partial  derivative  must  be  of  the  same 
sign  for  al_l  values  of  the  other  known  parameters.  This  may  sometimes  be 
shown  by  manipulation  of  the  derivative  expression  to  show  it  is  always 
of  one  sign  due  to  its  form.   In  other  cases  a  worst  case  technique  will 
be  sufficient  to  prove  a  constant  sign  of  the  derivative.   In  this  pro- 
cedure, if  the  derivative  is  assumed  to  be  always  positive  all  known  para- 
meters appearing  in  the  expression  which  tend  to  increase  the  derivative 
are  set  at  their  minimum  and  all  those  which  tend  to  decrease  the  deriva- 
tive are  set  at  their  maximum.   If  the  derivative  remains  positive  for  this 
worst  case,  it  can  be  assumed  always  positive.  This  procedure  has  the 
weakness  that  it  may  not  be  evident  whether  a  given  parameter  tends  to 
increase  cr  decrease  the  derivative.  Also,  in  general,  if  the  range  of 
parameter  variation  is  very  large,  the  derivatives  will  not  pass  this 
worst  case  test  and  the  parameter  must  be  divided  into  increments  for  in- 
vestigation. 

When  the  ranges  of  all  parameters  have  been  divided  into  increments 
to  be  tested  whether  numerous  increments  or  only  one  (end  points  type  of 
parameter)  the  computer  program  is  written.  This  program  solves  for  the 
output  using  each  increment  of  each  design  parameter  in  combination  with 


aach  increment  of  each  knows  parameter.  The  output  resulting  from  each 
set  of  parameters  is  tested  and  those  sets  which  result  in  output  which 
meets  specifications  are  stored  and  printed. 

The  output  of  the  computer  can  most  easily  be  interpreted  if  pre- 
sented in  graphical  form.   If  there  is  just  one  design  parameter  this 
is  not  necessary.  A  roost  useful  presentation  for  two  design  parameters 
is  a  plot  of  one  versus  the  other.   Starting  with  the  data  for  one  value 
of  a  selected  known  parameter;  for  each  value  of  one  design  parameter^,  the 
maximum  and  minimum  acceptable  values  of  the  other  are  plotted.  Connect- 
ing all  such  points  will  yield  an  area.  The  values  of  design  parameters 
corresponding  to  any  point  within  this  area  when  used  with  the  previously 
selected  known  parameter  will  result  in  a  circuit  which  will  meet  design 
specifications. 

This  procedure  is  repeated  until  area  plots  are  obtained  for  every 
increment  of  each  of  the  known  parameters.   In  each  case  one  design  para- 
meter is  plotted  versus  the  other  utilizing  the  same  scales.   If  all  these 
plots  are  overlayed,  the  area  which  is  common  to  the  acceptable  area  of 
all  the  graphs  is  the  final  design  area  desired.   If  a  point  in  this  final 
common  area  is  selected  and  the  corresponding  values  of  design  parameters 
are  used  in  a  circuit,  the  circuit  will  meet  specifications  for  all  known 
parameters. 

The  requirement  for  obtaining  an  area  graph  for  each  increment  of 
each  known  parameter  demands  much  tedious  plotting.  This  is  essentially 
the  procedure;  however,  it  is  possible  to  write  the  computer  program  so 
that  the  data  will  be  partially  processed  by  the  computer  before  printing 
thus  eliminating  some  of  the  plotting  and  overlaying.  This  programming 
method  is  most  easily  described  by  use  of  a  specific  example  and  is  de- 
monstrated in  a  later  chapter. 

7 


It  Is  quite  possible  that  when  the  acceptable  area  graphs  for  the 
increments  of  known  parameters  are  overlayed  there  will  be  no  acceptable 
area  which  is  common  to  all  the  graphs .  Thus  there  are  no  combinations 
of  design  parameters  which  will  guarantee  satisfactory  circuits  for  ran- 
dom known  parameters.  This  condition  can  occur  because  the  specifications 
on  the  output  are  too  strict.  This  may  be  stated  another  way  by  saying 
that  the  range  of  one  or  more  of  the  known  parameters  is  too  large.  The 
difficulty  may  be  corrected  by  relaxing  the  specifications  on  the  output 
or  dividing  one  of  the  known  parameters  into  two  groups.   If  the  trouble- 
some known  parameter  is  divided  into  two  groups,  each  of  which  has  half 
the  range  of  variation  of  the  parameter,  there  may  then  be  common  areas 
of  acceptable  values  for  these  half  ranges  and  all  other  known  parameters. 
In  this  case  then  there  are  in  effect  two  design  areas  rather  than  one  and 
the  divided  parameter  must  be  tested  and  separated  into  two  groups  before 
one  can  use  this  procedure  to  design  a  circuit. 

If  the  area  is  large,  it  may  be  possible  to  tighten  the  output  speci- 
fications to  be  met  or  to  add  another  specification.  Doing  either  of 
these  will  tend  to  reduce  the  size  of  the  area,.   By  inspecting  the  size 
of  the  area  resulting  from  the  initial  output  specifications,  a  decision 
can  be  made  as  to  how  the  specifications  can  be  altered.  By  inspecting 
subsequent  results  and  making  further  adjustments  a  balance  between  specif- 
ications and  available  combinations  of  design  parameters  can  be  achieved. 

The  technique  described  results  in  an  area  of  acceptable  design 
values.   If  the  designer  picks  a  set  of  design  values  from  in  the  final 
area,  he  is  assured  that  the  circuit  •sill  meet  oafcrmt  specifications. 
The  question  still  remains,  "If  all  the  combinations  in  the  area  meet 
specifications  is  there  not  one  best  combination  or  some  part  of  the  area 
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which  yields  better  circuits  than  other  parts?"  The  techniques  des- 
cribed up  to  the  point  will  not  answer  this  question. 

One  method  of  determining  the  best  part  of  the  final  common  area  is 
to  select  one  of  the  preliminary  area  plots.  This  plot  corresponds  to 
one  set  of  known  parameters.  Using  this  set  of  known  parameters  and  re- 
presentative values  of  design  parameter  combinations  throughout  the  area, 
one  calculates  the  value  of  one  of  the  output  quantities.  The  values  cal- 
culated are  plotted  at  the  corresponding  point  in  the  area.   Inspection  of 
the  plot  may  reveal  a  trend  in  the  value  of  the  output  parameter  and  in- 
dicate a  portion  of  the  area  where  this  value  best  approximates  its  opti- 
mum value.  This  process  may  be  repeated  with  other  sets  of  known  para- 
meters and  the  corresponding  area  until  it  is  established  that  the  trend 
holds  in  enough  of  these  preliminary  areas  that  it  may  be  assumed  to  hold 
in  the  final  area  plot. 

In  summary  the  technique  which  has  been  described  in  this  section 
results  in  a  plot  of  one  design  parameter  versus  the  other.  An  area  of 
acceptable  combinations  of  the  two  design  parameters  can  be  drawn  on  this 
plot.   If  a  combination  is  selected  from  inside  this  area  and  a  circuit 
constructed  using  this  combination  and  any  combination  of  the  known  para- 
meters used  in  the  circuit,  the  resulting  circuit  will  meet  the  desired 
output  specifications. 


3.  Description  of  Neon-photoconductor  Circuit. 

The  circuit  chosen  to  illustrate  the  technique  described  briefly  in 
the  preceding  chapter  is  a  neon  tube  -  photoconductor  circuit  which  is 
used  as  a  basic  building  block  in  the  design  of  computer  logic  circuits. 
The  basic  circuit  is  shown  in  figure  1. 


. 
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Figure  1. 

The  circuit  consists  of  a  neon  tube,  two  fixed  resistances,  Rl  and 
R2  and  a  photoconductor,  Rpl,  supplied  from  a  300  volt  source. 

In  the  circuit,  the  output  quantity  of  interest  is  the  light  of  neon 
tube  two.  The  input  is  light  from  a  similar  external  neon  tube,  called 
number  one.   With  no  input,  there  is  no  output;  the  very  high  dark  resist- 
ance of  the  photoconductor  which  is  in  the  mega  ohm  region  limits  current 
flow  and  keeps  the  voltage  drop  across  neon  number  two  low  enough  to  pre- 
vent its  firing.  With  input  light  from  neon  number  one,  the  resistance 
of  the  photoconductor  drops  radically  to  its  light  resistance  value  of  a 


few  ktlohms  and  neon  number  two  fires,  producing  light  output. 

The  light  output  from  neon  number  two  is  used  as  input  to  sub- 
se  -jcnt  circuits  so  it  is  essential  that  it  be  as  nearly  constant  as 
possible  for  stable  operation.  This  constant  current  characteristic 
is  not  possible  with  the  circuit  as  shown  since  the  resistance  of  the 
photoconductor  is  a  function  of  temperature.  As  the  circuit  operates, 
the  ambient  temperature  rises,  causing  the  photoconductor  resistance  to 
rise  and  the  neon  current  to  decrease.  A  proposed  method  of  compensation 
intended  to  minimize  this  temperature  effect  is  to  insert  a  similar  photu- 
conductor  in  the  shunt  branch  of  the  circuit  as  shown  in  figure  2. 
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Figure  2. 
In  this  compensated  circuit  the  output  current  is  more  nearly  con- 
stant despite  temperature  increase  because  of  the  feedback  nature  of  the 
output  light  shining  on  photoconductor  Rp2.  Now,  as  the  output  current 
tends  to  decrease  as  a  result  of  temperature  variation  changing  the  re- 
sistance of  Rpl,  the  light  output  decreases,  increasing  the  resistance  of 
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Rp2  and  shunting  more  of  the  available  current  through  the  neon  tending 
to  maintain  the  current  constants 

The  problem  then  is  to  design  a  circuit,  using  available  components, 
which  will  give  a  nearly  constant  output  light  or  current ,  throughout 
the  operating  temperature  range  of  25  to  55  degrees  Centigrade.  The 
supply  voltage  is  assumed  to  be  a  constant  300  volts.  The  neon  to  be 
used  is  a  General  Electric  NE2  type  which  has  a  firing  voltage  of  135  volts 
and  a  maintaining  voltage  of  from  60  to  100  volts. 

The  characteristics  of  the  photoconductors  employed  in  the  circuit 
are  more  complex  since  the  value  of  photoconductor  resistance  is  a  func- 
tion of  both  the  ambient  temperature  and  the  incident  light.  The  latter 
variation  is  more  easily  investigated  if  neon  current  rather  than  light 
is  considered  the  determining  factor.  Mathematical  relationships  governing 
photoconductor  resistance  variation  with  temperature  and  with  current  have 
been  determined  from  experimental  data. 

The  temperature  dependence  of  resistance  can  be  obtained  by  consider- 
ing figure  3  where  log.n  of  resistance  in  kilohna   is  plotted  versus  tempera- 
ture in  degrees  Centigrade  for  three  typical  photo  conductors.  The  curves 
are  all  for  a  constant  incident  light  corresponding  to  two  milliamperes  of 
neon  current.   Since  the  data  plots  as  a  straight  line  on  semi- log  paper 
the  variation  is  of  the  form  Kp~  r\  <  if  curve  II  is  selected  to 

determine  the  constant, 

PPG*0    -    A  <?***   =    J.3  a/**  (i) 
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The  coefficient  A  and  the  mathematical  relationship  expressing  re- 
sistance variation  with  neon  current  can  be  obtained  by  considering 
figure  4.   In  this  curve,  photoconductor  resistance  in  kilohms  is  plot- 
ted versus  current  in  the  exciting  neon  expressed  in  milliampares;  the 
curves  are  all  for  a   constant  temperature  of  40°  C.   The  three  curves 
A,  B,  and  C  are  for  tbree  representative  photoconductors,  A  and  C  show 
the  extremes  of  variation  and  curve  B  is  considered  typical.   In  attempt- 
ing to  express  the  photoconductor  resistance  as  a  function  of  light,  it 
was  noted  that  curves  A,  B  and  C  are  all  of  the  same  mathematical  form 
and  are  related  by  a  constant.  This  magnitude  multiplying  factor  will 
be  defined  as  Z.   This  relationship  is  best  illustrated  by  table  I,  figure 
5. 

The  current-resistance  characteristic  of  any  photoconductor  can  be 
obtained  by  deriving  a  mathematical  expression  which  will  fit  one  of  the 
curves  giving  the  proper  form  of  variation  and  then  multiplying  the  en- 
tire expression  by  a  constant  to  fix  its  magnitude. 

Curve  B,  selected  as  typical  will  be  the  one  used  to  determine  the 
mathematical  form  of  variation.   If  this  expression  is  assumed  to  have  a 
magnitude  multiplying  factor  of  unity,  reference  to  figure  4  shows  that 
in  order  to  cover  the  range  from  curve  A  to  curve  C  the  magnitude  multiply- 
tng  factor  must  have  values  of  from  .4,  for  curve  A  to  2.5  for  curve  C  . 

The  first  approximation  to  curve  B  which  was  investigated  is  a 
hyperbolic  function  where  K1  is  a  constant  to  be  determined. 

Since  a  neon  current  of  2  ma.  will  be  used  as  the  mean  point  about  which 

to  design,  the  hyperbola  constant  was  determined  to  fit  curve  B  exactly 

3 
at  that  point.   From  tie  curve  at  I  ■  2.0  ma.,  Rp  ■  4.388  x  10  ohms. 
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Tabl 

.     i 

Current 
in  ma. 

Rp 

Curve  A 
in  Kohm 

Rp 

Curve  B 
in   Kohm 

Rp 
Curve  C 

in   Kohm 

Rp(A) 

Rp    (B) 

Rp(ci 

Rp    (B) 

1.00 

A.  48 

12.66 

30.8 

.354 

2.43 

1.25 

3.38 

9.97 

22.0 

.370 

2.43 

1.50 

2.78 

6.71 

17.0 

.414 

2.53 

1.75 

2.06 

5.32 

13.86 

.388 

2.60 

2.00 

1.72 

4.388 

11.22 

.394 

2.54 

2.25 

1.46 

3.79 

9.28 

.396 

2.44 

2.50 

1.26 

3.37 

7.88 

.374 

2.34 

2.75 

1.11 

3.04 

6.76 

.363 

2.22 

3.00 

1.00 

2.89 

5.82 

.346 

2.02 

Figure  5 

a    *    _£..  lit 

.  .     ^  p 
Incorporating  the  effect  of  temperature, 

3 
Since  the  value  of  P.  s  4.388  x  10  ohms  from  curve  B  corresponds  to 

I  *  2ma.  and  T  ■  40°C  the  constant  8.776  Kj  can  be  determined 

?<776      /So.      -   -     ? 
The  hyperbolic  expression  with  temperature  effect  included  approximat  ing 

curve  B  is  then 

X 

To  include  the  fact  that  the  magnitude  of  the  resistance  can  vary 
from  curve  A  values  to  curve  C  values,  the  magnitude  multiplying  factor 
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Z  must  be  included.  The  complete  resistance  expression  for  the  photo- 
conductor  is  then 

n   ■ 

In  the  circuit  being  analyzed,  Rpl  is  supplied  with  exciting  light 
from  an  external  neon,  neon  number  1,  which  is  assumed  to  provide  a  con- 
stant light  intensity  corresponding  to  a  neon  current  of  2  ma.  There- 
fore •         O  ^ 

The  compensating  photo  conductor  Rp2  is  excited  by  light  fed  back 
from  the  output  neon,  neon  number  2.   This  light  corresponds  to  the  out- 
put current  I-.  Therefore : 
V  2  C/S 

£7      <?   =    -Y»J?£  &      <?' 

The  plot  of  the  hyperbolic  approximation  to  curve  B  shown  in  Figure 
6  indicates  that  the  hyperbolic  curve  fits  experimental  data  exactly 
only  at  I  ■  2,0  ma.   Since  circuit  performance  with  currents  from  1.5 
ma.  to  2.5  ma.  is  of  interest,  the  use  of  the  hyperbolic  approximation 
introduces  some  error  especially  at  the  lowest  value  of  1.5  ma.  where 
the  error  is  12.8%. 

The  expression  for  photoconductor  resistance  is  used  in  subsequent 
calculations  and  the  simplicity  of  the  hyperbolic  expression  is  an  ad* 
vantage  when  used  in  more  complex  expressions.   However,  a  more  accurate 
curve  fit  over  the  entire  1.5  to  2.5  ma.  current  range  was  desired. 

The  second  approximation  to  curve  B  which  was  investigated  was  de- 
rived by  assuming  the  curve  to  be  of  the  forms 

A  PrlX  +    8X     f    C  £P    f-P  --  o 

where  A,B,C,  and  D  are  constants  to  be  determined.   Appendix  I  shows  the 
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calculations  leading  to  the  final  form 


from  which 


Including  the  effect  of  temperature  as  in  the  hyperbolic  approxi- 


mation: /  ,  \ 

1 


In  the  circuit  being  analyzed,  Rl  is  supplied  by  light  from  neon  one 
which  has  a  constant  current  of  2.0  ma.   therefoia? 


■ 


resistance  have  such  large  ranges  of  variation,  an  output  current  of  from 
1.5  to  2.5  ma.  will  meet  specifications.  Rl  and  R2  must  be  chosen  such 
that  this  specification  is  met  with  any  neon  tube  and  any  photoconductor 
having  characteristics  within  the  ranges  discussed  above.  The  only  assump- 
tion being  that  both  photoconductors  have  the  same  value  of  magnitude  multi- 
plying factor,  Z.  A  further  specification  is  that  for  any  finished  circuit 
which  will  have  one  value  of  neon  voltage,  V,  and  one  value  of  Z,  the  vari- 
ation in  current  due  to  temperature  change  from  25°C  to  55°C  shall  not  be 
greater  than  .15  ma„  The  effect  of  imposing  stronger  restrictions  on  the 
output  will  be  investigated  and  the  results  compared. 

The  upper  limit  on  the  values  of  Rl  to  be  investigated  can  be  estab- 
lished by  setting  the  current  II  »  1.5  and  the  voltage  V  *  60, 
then 

Rl    cannot  be  as  conveniently  set  since  its  value  depends  upon  the  value 

of  R2  which  has  not  yet  been  determined.  A  rough  initial  estimate  can  be 

obtained  by  assuming  no  current  flow  through  the  shunt  branch  and  ignoring 

the  contribution  of  Rpl  in  comparison  with  Rl  then 

<P«?  X  /c 
Jv. 

In  view  of  the  approximations  made  here,  a  more  pessimistic  value  of  20 

kilohms  should  probably  be  chosen.  Experience  has  shown  that  a  range  of 

3  3 

from  40  x  10  to  130  x  10  ohms  for  Rl  will  adequately  cover  the  possible 

values  of  Rl  to  have  the  circuit  meet  output  specifications. 

The  maximum  value  of  R2  is  indeterminant  since  the  effect  of  increas- 
ing R2  is  to  shunt  more  current  through  the  neon  tube.   It  can  be  seen  that 
if  II  is  2.5  ma.  or  less,  R2  can  be  increased  without  bound  and  the  current 
12  will  never  exceed  the  output  specification.  The  minimum  value  of  R2  is 
also  indeterminate.  In  the  uncompensated  circuit  a  minimum  value  of  R2  can 

20 


be  determined  because  a  lower  limit  exists  below  which  the  neon  voltage 

will  not  be  above  the  135  volts  necessary  for  firing.   In  the  compensated 

circuit  however,,  due  to  the  very  high  dark  resistance  of  the  compensation 

photoconduceor,  the  neon  will  fire  for  any  value  of  R2.   The  magnitude 

of  Rl    .      was  selected  as  an  initial   lower  limit   for  R2  also.     An  initial 
min 

upper  limit  of  700  K  ohms  was  selected.   Subsequent  experience  showed 

3  3 

that  a  range  of  R2  from  40  x  10  to  530  x  10  ohms  was  sufficient  to   yield 

the  desired  design  curves. 
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4.   Preliminary  Ci*      .raalysLs. 

This  section  deals  with  the  application  of  the  design  technique  deve- 
loped in  Section  2  to  the  aeon-photoconductor  circuit  which  has  been  pre- 
viously described, 


Figure  7. 
The  output  quantity  of  interest  is  the  light  from  neon  number  2. 
As  previously  noted,  the  output  is  mora  easily  expressed  as  the  current 
through  the  aeon,  an  equally  good  measure  of  performance.  The  technique 
first  requires  an  equation  expressing  the  output  as  a  function  of  the 
other  parameters. 


£1   +   £\ 
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0      \ 
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If  the  hyperbolic  approximation  is  used  for  photoeoaduetor  resistance, 

2  f 
This  equation  eats  ha  easily  manipulated  into  the  form  of  a  quadratic 

where 


0*j*MXy  T(v-  *^C'f3rJ 

-1 

Then  by  the  quadratic  equations 

Inspection  of  the  coefficients  A,B,  and  C  shows  that  since  A  is 

always  positive  and  C  is  always  negative  (the  maximum  value  ©f  V  is  100) 
the  quantity  resulting  from  evaluating  the  radical  will  he  greater  than 
B.  Since  current  12  cannot  physically  he  in  the  direction  opposite  to 
that  assuwedj  a  negative  answer  for  12  has  no  significance  and  the  equAT 
tion  to  be  solved  with  the  computer  is 

If  the  second  order  approximation  for  photoconductor  resistance  is 


used: 


2l  +  AllX/o'*£'c"2'' 
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This  squat  ion  can  also  be  manipulated  into  the  quadratic  form; 
where  the  coefficients  are  somewhat  more  complex, 


4       _  r      >Cll 


+    iHZ  X  /      ,  01' 

1, 

■    -  emC,i3T  |        */cj  v     -, 


B^     j*;*y.*x 


+  /-Hi  x/.       //ei 


The  same  quadratic  formula  must  then  be  solved: 


In  order  to  determine  whether  any  of  the  known  para/nefe«? rs  fall  in- 
to the  category  where  only  the  maximum  and  minimum  values  need  to  be  in- 
vestigated, partial  derivatives  of  12  with  respect  to  V,  Z,  and  T  must 
be  taken  and  checked  to  see  if  they  are  of  constant  sign. 
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Inspection  shows  that  the  quantities  Aj,  " — ~  '  and  — —  are  always 
positive,  C  and-        always  negative  and  the  sign  of  S  may  be  plus 
or  minus  depending  upon  the  magnitude  of  the  negative  term 
compared  with  the      ive  term 

The  large  number  of  terms  in  the        expression  make  it  unwieldly 
and  no  nsnipulation  which  would  establish  tha  sign  of        suggests 
itself.  Therefore,  the  second  approach  must  be  used.  The  slope  will  be 
assumed  negative  and  the  "worst  case"  values  of  Rl,  R2,  T,  Z,  and  V  which 
will  tend  to  make  the  slope  positive  will  be  inserted.  The  expression  with 
these  values  inserted  is: 

ii" 


e 


'•j  v        =  \  ~3  \ 


t~ 


■ 


The  expression  evaluated  is  a  +  .0374.  A  positive  answer  shows  that 
the  initial  assumption  of  a  negative  slope  may  be  incorrect.   The  informa- 
tion gained  from  this  test  is  not  conclusive,  since  it  is  improbable  that 

25 


the  worst      values  of  all  the  parameters  will  occur  simultaneously, 
However s  since  a  positive  slope  was  calculated,  the  pessimistic  alter- 
native must  be  *       id  the  temperature  must  be  investigated  in  in- 
crements. A  similar  worst  case  analysis  on  <*  ■*"*'     leads  to  the  same 
conclusion, 

The  variation  of  12  with  V  will  be  investigated  next. 

A   -J  C 


J?.  J  i  X  /        &  '  +•  J?  1  1 


« — ^  *  tA  s*„ 


°>  v         j/ll/  t-  B 


J  f 5      w  • 
Since  — — ^  and      are  always  greater  than  zero  ■     will  be  nega- 

tive  as  long  as 
Since 
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an-.        this  last  Inequality  holds  and  thus  the  in- 
itlal  assumption  of  a  negatl         is  txme.   Since        is  always 
negative,  only  the  maximum  and  minimum  values  of  V  are  investigated. 

When  a  similar  analysis  of  these  partial  derivatives  is  made  using 
the  more  accurate  second  order  approximation,  it  is  found  that  ^     , 

and  all  have  possible  changes  in  sign  and 

therefore  when  this  approximation  is  used  all  three  parameters  must  be 
investigated  using  increments.   This  condition  is  then  one  disadvantage 
to  using  the  second  order  approximation. 

When  it  has  been  determined  that  a  parameter  must  be  investigated 
using  increments,  the  magnitude  of  the  increments  must  be  chosen.  As 
pointed  out  in  Section  2,  the  number  of  increments  selected  is  a  com- 
promise between  accuracy  and  keeping  the  number  of  calculations  reason- 
able. The  initial  choice  of  increments  was  as  follows: 
T  =  25,  35,  45,  55, 
Z  *  .4,  ,9,  1.4,  1.9,  2.5, 

V  =  60,  100,  for  hyperbolic  approximation, 

V  =  60,  70,  80,  90,  100,  for  second  order  approximation, 
Rl  »  in  10  kilohm  increments, 

R2  =  in  20  kilohm  increments. 

Plots  made  using  these  increments  showed  that  parts  of  the  range  of 
Rl  caused  the  value  of  R2  to  increase  rapidly  with  possible  inaccuracy  in 
the  plot.  The  reason  for  this  rapid  increase  is  explained  in  Section  6. 
To  eliminate  this  possible  inaccuracy  the  increments  of  Rl  were  altered 
so  that  in  the  questionable  region  the  increments  were  taken  much  smaller 
(2  Kilohm  increments). 

The  problem  has  at  this  point  been  reduced  to  the  solution  of  a 
quadratic  equation  with  coefficients  made  up  of  the  five  parameters  Rl, 
R2,  T,  Z  and  V.   The  equation  is  to  be  solved  for  all  combinations  of 
each  of  the  increments  of  the  five  parameters. 
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3.   Digital  Computer  Program  Technique 

The  solution  of  a  quadratic  equation  whose  coefficients  are  func- 
tions of  five  variables  and  where  the  number  of  solutions  is  the  product 
of  the  parameter  increments  is  ideally  suited  to  a  digital  computer.  A 
parameter  increment  is  defined  as  the  number  of  values  that  a  variable 
will  be  assigned  spanning  its  entire  range  of  variation.  The  number  of 
combinations  of  parameters  to  be  examined  in  this  case  will  be  denoted  by 

N  where: 

(  0 1   \  y /?*/  )  X  /  "7l    H  c 

Programming  of  a  digital  computer  in  machine  language  requires  an 
elaborate  and  complex  sequence  of  instructions  to  specify  every  step  of 
an  operation.   With  the  advent  of  the  FORTRAN  processor  in  late  1956, 
the  intricacies  of  programming  were  no  longer  a  stumbling  block.  Al- 
though the  FORTRAN  system  was  developed  by  J.  W.  Backus  at  IBM,  it  has 
found  universal  acceptance  for  all  types  of  digital  computers.   The  FORTRAN 
processor  has  made  it  possible  for  scientists  and  engineers  to  apply  the 
language  of  ordinary  mathematics  to  their  individual  problems. 

Quoting  from  a  paper  prepared  by  Backus  and  his  associates,  "The 
purpose  of  FORTRAN  was  to  reduce  by  a  large  factor  the  task  of  preparing 
scientific  problems  for  IBM's  next  large  computer,  the  704.   If  it  were 
possible  for  the  704  to  code  problems  for  itself  and  produce  as  good  pro- 
grams as  human  coders  (but  without  the  errors),  it  was  clear  that  large 
benefits  could  be  achieved.   For  it  was  known  that  about  2/3  of  the  cost 
of  solving  most  scientific  and  engineering  problems  on  large  computers 
was  that  of  problem  preparation.   Furthermore,  more  than  90  per  cent  of 
the  elapsed  time  for  a  problem  was  usually  devoted  to  planning,  writing 
and  debugging  the  program.   In  many  cases  the  development  of  a  general 
plan  for  solving  a  problem  was  a  small  job  in  comparison  to  the  task  of 
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devising  and  c<>       tiiite  procedures  to  carry  out  the  plan.  The  goal 
of  the  FORTRAN  project  was  to  enable  the  programmer  to  specify  a  num- 
erical procedure  using  a  concise  language  like  that  of  mathematics  and 
obtain  automatically  from  this  specification  an  efficient  704  prog 
to  carry  out  the  procedure.   It  was  expected  that  such  a  system  would 
reduce  the  coding  and  debugging  task  to  less  than  one- fifth  of  the  job 
it  had  been." 

FORTRAN  (derived  from  formula  translation)  is  a  language  that  great- 
ly resembles  the  language  of  ordinary  mathematics.  Symbols  such  as,  **, 
*»  /»  +»  "»  ifl  FORTRAN  are  digested  as  exponentiation,  multiplication, 
division,  addition,  and  subtraction.  Mathematical  operations  such  as 
square  roots,  logarithms,  trigonometric  and  hyperbolic  functions  are 
available  in  the  computers  library  of  routines  and  used  merely  by  re- 
ferencing the  appropriate  operation  by  abbreviated  name. 

The  Control  Data  Corporation  1604  digital  computer  was  married  to 
the  FORTRAN  processor  in  late  1961.  The  results  of  the  merger  produced 
a  high  speed  computer  that  could  be  commanded  with  relative  e^se.   Suc= 
cessful  results  with  FORTRAN  programming  can  best  be  achieved  by  follow* 
ing  these  basic  rules  in  the  order  shown. 

1.  Block  diagram  -  (flow  chart) 

2.  Mathematical  and  logical  operations  to  be  performed 

3.  Transformation  to  FORTRAN  language. 

An  example  of  a  basic  flow  chart  for  the  solution  of  the  quadratic 
equation  using  the  hyperbola  approximation  is  illustrated  in  figure  8. 
The  idea  behind  this  program  is  to  see  the  effect  on  the  neon  current  as 
the  temperature  changes  for  any  one  set  of  known  and  design  parameters. 
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SET 

V    -  V    ♦   1 


5ET 

Z   .  I   .   1 


jET 
HI    -  HI    t   1 


SLT 

t  -  r  +  i 


EMER 
A-    IN   UT   J.   LUT-UT    FGHMATS 

B-   PAHA'ETER   VARIABLES   -    INPUT   UATA 


SET   VOLTAGE    T(l    IMTIAL   VALUE 


SLT    I    TCJ    IMTIAL   VALUE 


SET   HI    TO    IMTIAL   VALUE 


d£T   n2   TO    IMTIAL   VALUE 


SET  CMAX  \   emu;  TO  OUTSIUE  VALUES 


SET  T  TO  IMTIAL  VALUE 


EALCULATE  NEON  CLHnEM 


IS  EIHEUIT  AEEEPTAQLE? 


STORE  MAX  L   MIN  CURRENT 


HAVE  ALL  T  INCREMENTS  UEEN  THIE07 


YES 


OELTA  »  MAX  -  MIN 


(SITE  OUTPUT  MAGNETIC  TAPE 
V  Z  HI  H2  DELTA 


HAVE   ALL   H2    INCREMENTS   OEEN   TRIED? 


HO 


*LL   Hi    [NCREI  ENTS   BEEN   TRIED? 


N( 


YES 


IAVE   ALL   7.    INCREMENTS  OEi.fi   THIEO? 


YES 


HAVE  ALL  V  INCREMENTS  SEEN  TRIED? 


NG 


PROGRAM  l'.j 


Fig.    8 
3  0 


SET 
R2  -  R2  +  1 


The  fiow  chart  indicates  the  sequence  of  operations  to  be  performed 
with  alternate  paths  dependent  on  logical  tests*   Each  block  in  the  flow 
chart  may  represent  a  few  calculations  or  logical  interrogations.  Es- 
sentially the  flow  chart  in  figure  8  systematically  takes  every  possible 
combination  of  circuit  parameters  and  tests  the  output  current  to  see  if 
it  meets  specifications.   If  a  failure  occurs  the  design  parameter  is 
altered  and  the  calculations  and  logical  tests  are  started  over  again. 
When  the  entire  range  of  temperature  variation  has  been  investigated  and 
yielded  an  acceptable  circuit,  the  design  and  known  parameter  set  which 
was  responsible  for  the  acceptable  circuit  for  the  entire  temperature  varia- 
tion is  recorded  on  magnetic  tape.   In  addition  to  the  parameter  set  the 
variation  in  neon  current  is  also  placed  on  the  magnetic  tape.  This  pro- 
cess continues  with  a  new  parameter  combination  until  all  parameter  combin- 
ations have  been  investigated. 

The  mathematical  calculations  to  be  performed  in  the  solution  of  a 
quadratic  equation  consist  of  evaluating  the  coefficients  which  are  vari- 
ables dependent  upon  the  parameter  set  under  consideration  and  then  using 
the  quadratic  formula  to  evaluate  the  neon  current.   Since  in  the  present 
problem  the  argument  of  the  square  root  is  always  positive  and  greater 
than  the  coefficient  B,  only  the  positive  root  need  be  calculated.   Logi- 
cal tests  to  be  performed  consist  of  determining  acceptable  circuit  designs 
by  employing  the  given  circuit  specifications  as  test  criteria.  Other  logi- 
cal tests  are  used  in  determining  maximum  -  minimum  values  and  program  loop 
completions. 

The  transformation  from  mathematical  to  FORTRAN  statements  is  accom- 
plished next.   Input  to  the  1604  digital  computer  can  be  by  either  punch- 
ed cards  or  magnetic  tape.  Utilization  of  punched  cards  is  superior  at 
this  time  as  cards  which  are  in  error  can  be  easily  replaced,  and  addi- 
tions or  deletions  made  as  necessary.  The  complete  program  in  FORTRAN 
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appears  as  in  figure  99  Each  line  in  the  program  represents  a  single 
punched  IBM  card.  If  this  same  program  had  been  done  in  machine  lan- 
guage, these  thirty  odd  FORTRAN  cards  would  be  replaced  by  about  five 
hundred  machine  language  cards. 

If  the  FORTRAN  program  when  placed  in  memory  is  complete  and  logi- 
cal, the  program  will  compile.   If  compilation  is  not  accomplished  an 
error  printout  will  result.  The  error  printout  categorizes  the  errors 
that  have  escaped  detection  by  the  programmer  and  lists  these  errors  by 
FORTRAN  statement  number.  When  the  error  is  recognized  altetations  to 
the  punched  cards  can  easily  be  accomplished  and  the  program  compilation 
retried  until  successful.  Writing  output  on  magnetic  tape  is  a  separate 
entity  from  the  computer  proper  and  occurs  simultaneously  with  computer 
calculations.  The  magnetic  tape  in  turn  is  read  out  on  an  off  line 
printer  so  as  to  conserve  computer  time.  An  example  of  the  digital  com- 
puter printout  as  obtained  from  the  off  line  printer  appears  as  in  figure 
10. 

Prior  to  proceeding  with  these  digital  computer  answers  a  mathemati- 
cal check  is  required  to  see  whether  the  problem  is  being  solved  correct- 
ly. A  desk  calculator  should  be  used  rather  than  a  slide  rule  to  obtain 
enough  significant  figures  for  comparison  purposes.  When  a  calculator 
check  has  been  made  for  one  set  of  parameters  and  dependability  of  the 
computer  program  proven,  then  the  accuracy  of  the  remaining  output  can  be 
relied  upon.  A  similar  FORTRAN  program  for  the  uncompensated  neon  photo- 
conductor  circuit  is  illustrated  in  figure  11.  This  program  is  needed 
for  a  comparison  of  its  results  with  those  of  the  compensated  circuit. 
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PROGRAM  BASIC         .  „ 

DIMENSION  V(3),R1(  10),'ft2(  17)  ,T(4),Z(4) 
IOO  F0RMAT13F4.0) 

READ  I00(V( I ). 1=1,3) 
200  FORMAT  I  10F7.0)         ... 

READ  200  (R1( J),  J=1, l6l 
300  FORMAT  (  11F7.0) 

READ  300(R2U),K»li  17) 
400  FORMAT  (4F3.0) 

READ  400  (T(L),  L=1,4) 
500  FORMAT  (4F3. 1 ) 

READ  500  (Z(M),M=1,4) 

DO  16  1=1,3 

DO  15  M=1,4 

DO  14  J  =  1,10 

DO  13  K=l,  17 

CMIN1=3. 

CMAX1=1.. 

DO  12  L  =  1,1* 

A  =  (R1( J)*R2(K)+2. 1 1E3»Z(M)»R2(K)»EXPF{ .0183*T(L) ) 1/1.E6 

E  =  2.  1  1E3»Z(M)*EXPF{.0183»T(L)  )»(V(  I  )  +  (2.»R  I  (  J)  )/I.E3KV(  I  l»Rl  I  J) 

F  =  8.904E3»Z(M)*Z(M)«EXPF{ .0366*T(L  )  )+R2(K)*(V( I  )-300.) 

8  =  (E+F1/1.E3 

C  =  4.22*Z(M)»EXPF( .0183»T(L) )»( V( D-300.) 

CTWO  =  <-B+SQRTF( (B«B)-(4.*A*C) ) )/(2.*A) 

IF  (CTWO-1.5) 13, 13,8 
.   8  IF  (CTW0-2.5)9, 13, 13 

9  CMAX1=MAX1F(CMAX1,CTW0) 

CMIN1=MIN1F(CMIN1,CTW0) 

12  CONTINUE 
DELTA=CMAX1-CMIN1 

600  FORMAT (F5.0.2F1 3.0, F8.  1.F13.5) 

WRITE  OUTPUT  TAPE  4 , 600 , V( I ) ,R 11  J ) ,R 2 ( k ) , Z ( M ) , DE LTA 

13  CONTINUE 
U  CONTINUE 

15  CONTINUE 

16  CONTINUE 
END  FILE  4 
PAUSE  1 
END 

END 

60.  80.100. 

40Q00.  50000.  60000.  70000.  80000.  90000.  1 00000. 1  10000. 120000. 1 3U000. 
2 00000. 220000. 240000. 2 60000. 2 80000.300 000.  320000  .  3U0  00  0.  36  0000  .  V 30000  .  4000  00 
420000. 440000. 460000. U 80000. 500000.520000. 
25.35.45.55. 
2.51.71.0  .4 


Fig.  9 
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PROGRAM  COMPA 

DIMENSION  V(  3),  R.M.I  0)  ,  R?  (  1  7  )  ,  T  (  4  )  ,  7(  4  ) 

F0RMATI3F4.0)     <  ' 

READ  100(V( I ), I=<1,3) 

FORMAT  (  IOF7.0) 

READ  200  (R  I  (J),  J= 1 ,  10) 

FORMAT  (  1  IFf.O) 

READ  300(R2<K),K=1,.17) 

FORMAT  (4F3.0) 

READ  400  (T(L  ),  L=;1  %M 

FORMAT  (4F3.  1  )        ' 

READ  500  (Z(M),M^1,4) 

DO  16  1*1,3 

DO  15  M=1f 4 

DO  14  J  =  1,10 

DO  13  K=l,  17 

CMIN 1=3. 

CMAX1=1. 

DO  12  L  =  1,4 
0CTWO=( (300.-VI  I  )  )/(Rl( J) +2.  I  1  E  3 »  Z ( M ) *EXPF( . 01H3*T(L) )  )-( VI  I  )/ 
1R2(K)))»1.E3 

IF  (CTW0-1  .5)  13,  13,0 

8  IF  (CTW0-2.SJ9, 13, 13 

9  CMAX  1=MAX1F(CMAX l,CTW0) 
CMIN1  =  MIN1F(CMIN  I.CTWO) 

12  CONTINUE 
DELTA=CMAX 1-CMIN1 

600  F0RMAT(F5.0,2F I  3. 0 , F8 . 1 , F 1 3 . 5 ) 

WRITE  OUTPUT  TAPE  4 ,  600,  V  (  I  )  ,  R  1  (  J  )  ,R  2  (  K  )  ,  Z  (  M  )  ,  DF  LTA 

13  CONTINUt 

14  CONTINUt 

15  CONTINUE 

16  CONTINUE 
END  FILF  4 
PAUSE  1 
END 

END 

60.  80.100. 

40000.  50000.  60000.  70000.  80000.  90000.100000.110000.120000.130000. 
200000.220000.240000.2  60000.28  0000.300  000. 32  0000.314  000  0.360  000.380  000.40  0000, 
420000.440000.460000.4  80000.500000.520000. 
2.51.71.0  .4 


Fig.  11 


The  main  con1'  il   program  illustrates  one  of  the  prime  advantages 
of  the  digital  computer  in  evaluating  circuits  which  have  known  para- 
meters of  the  type  requiring  incremental  analysis.  The  flow  chart  of 
the  main  program  called  neon  or  light  which  uses  the  hyperbola  or  second 
order  approximation  is  illustrated  in  figure  12.  Here  again  every  possi- 
ble combination  of  circuit  parameters  are  systematically  tested  to  deter- 
mine if  they  meet  specifications.  With  each  set  of  design  parameters 
the  known  parameters  requiring  increments  are  combined  in  every  possible 
manner.   If  every  combination  of  these  incremental  known  parameters  pass- 
es all  logic  tests  with  one  set  of  design  parameters,  then  this  combina- 
tion of  design  parameters  is  recorded  on  magnetic  tape  along  with  the 
maximum  and  minimum  neon  current  for  the  particular  voltage.   This  signi- 
fies that  the  present  design  parameters  being  considered  have  produced  a 
circuit  which  has  successfully  passed  all  of  the  logic  tests  for  the  un- 
predictable known  parameters  Z  and  T.  The  FORTRAN  programs  for  neon  and 
light  are  depicted  in  figures  13  and  14. 

Each  set  of  output  data  in  figure  15  is  a  successful  design  combina- 
tion.  In  other  words  for  a  neon  tube  maintaining  voltage  of  60  volts, 
a  series  resistance  of  90  K  ohms,  and  a  shunt  resistance  of  200  K  ohms, 
the  neon  current  will  only  vary  between  2.32944  and  2.01816  miHiamperes 
for  the  hyperbola  approximation  and  between  2.32938  and  2.01802  millamper- 
es  for  the  second  order  approximation.  This  variation  in  neon  current  in- 
cludes any  variation  in  the  characteristics  of  the  photoconductor  due  to 
temperature  or  light  intensity  within  the  assumed  possible  limits.  The 
last  two  columns  of  the  output  data,  maximum  and  minimum  neon  currents  for 
the  chosen  design  parameters,  have  no  immediate  use  in  the  synthesis  of 
the  circuit.   Since  these  currents  had  to  be  calculated  to  evaluate  the 
success  of  the  circuit,  it  was  an  easy  task  to  store  them  in  memory  and 
print  out  the  maximum  and  minimum  values  for  each  design  parameter  set. 
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5ET 

V  •  V  +  1 


5ET 

Rl  =  Rl  +  1 


5ET 
Z  =  Z  +  1 


ENTER 
A  -  INPUT  1  DUTPUT  FORMATS 
0  -  PARAMETER  VARIADLE'j  -  INCUT  DATA 


dET   VLILTAGE    TO    INITIAl     i/ALUI 


SET   HI    TO    INITIAL   VALUE 


SET   R2   TO    INITIAL   VALUE 


SET   CHAX   S.   CMIN   TO   OUTSIDE   VALUES 


SET   Z    TO    INITIAL   VALUE 


CALEULATE   OUTPUT   CURRENT    (CT1) 
FOR   1   3t   TEMPERATURE 


IS   CIRCUIT   ACCEPTABLE? 


YES 


NO 


CALCULATE  OUTPUT  CURRENT  (CTN) 
FOR  N  th  TEMPERATURE 


IS  CIRCUIT  ACCEPTABLE? 


YES 


NO 


STORE  MAXCT  &  MINCT 


IS  (MAXCT  -  MINCT)  ACCEPTADLE? 


YE'.: 


NO 


RETAIN  MAX  OF  MAXS  4.  MIN  OF  MINS 


HAVE  ALL  Z  INCH  UEEN  TRIED? 


NO 


YES 


WHITE  OUTPUT  MAGNETIC  TAPE 
V,  Rl,  R2,  CMAX,  CMIN 


HAVE  ALL  R2  INCR  OEEN  TRIED? 


NO 


HAVE  ALL  Rl  INCR  BEEN  TRIED? 


X 


HAVE  ALL  V  INCR  OEEN  THIEO? 
1  NO 


YES 


PROGRAM  END 


SET 

R2  =  H2  +  1 


Fig.    12 


PROGRAM  NEON 

DIMENSION    VI3),RII2I ),R2(34) ,Z(5) 
100    FORMAT13F4.0) 

READ  100(V( I ), 1=1,3) 
200  FORMAT  (  1  1F7.0) 

READ  200  (R  I  (J) , J=l ,21 ) 

READ  200  (R2(K),K=1 ,34) 
500  FORMAT  (5F3. 1 ) 

READ  500  (Z(M),M=1f5) 

DO  16  1=1,3 

DO  15  J=1,21 

DO  14  K=1,34 

CMIN 1=3. 

CMAX1=1. 

DO  13  M=l,5 

A1=(R1(J)#R2(K)+2.11E3*Z(M)*R2(K)*EXPF( .4575) )/1.E6 

E1=2.11E3*Z(M)«EXPF(.45  75)»(V(I)+(2.*R1(J))/1.E3)+V(I)«R1(J) 

F1=8.904E3*Z(M)»Z(M)«EXPF(  .  9  1  50  ) +R2( K ) * ( V(  I  )-300.) 

B1=(E1+F1 )/1.E3 

C1=4.22»Z(M)*EXPF( .4575) *( VI  I )-300.) 

CT1=(-B1+SQRTF( (B1*B1)-(4.*A1*C1)))/(2.*A1) 

IF  t CT  1-  1  .5  )  14,4,4 

4  IF(CT1-2.5)   5.5, 14 

5  A2=(R1  ( J)*R2(K )+2.  1  1E3*Z ( M  )  *R2 ( K ) *EXPF ( .6405) )/1.E6 

E2=2. 11E3»Z(M)*EXPF( .6405)*( V( I)  +  (2.»R1(J))/1.E3)+VU)*R1(J) 

F2=8.904E3*Z(M)*Z(M)*EXPF(  1.281  )+R2(K)»(V( I  )-300. ) 

B2=(E2+F2)/ 1.E3 

C2=4.22*Z(M)*EXPF( ,6405)*(V( IJ-300.) 

CT2=(-B2+SQRTF( ( B2*B2 )- ( 4. »A2*C2 ) ) ) / ( 2. «A2 ) 

IF(CT2-I.5)  14,6,6 

6  IFICT2-2.5)   7.7,14 

7  A3=(R1( J)*R2(K)+2. 1  I E3*Z ( M ) »R2 ( K ) *EXPF ( . 8235 ) ) / 1 . E6 

E3  =  2. 1 1E3*Z(M)*EXPF( .8235  )«{ V(I)  +  (2:*R1(J))/1.E3)+V(I)*R1(J) 

F3=8.904E3«Z(M)»Z(M)*EXPF( 1 . 64  7 ) +R2( K ) # ( V( I )-300. ) 

B3=(E3+F3)/1.E3 

C3  =  4.22*Z(M)»EXPF( .8235 )•{ VI 1 1-300. ) 

CT3=(-B3+SQRTF( ( B3*B3 ) - ( 4. » A3*C3 ) ) ) / ( 2. »A3 ) 

IF(CT3-1.5)  14.8,8 

8  IFUT3-2.5)  9,9,  14 

9  A4=(R1( J  )»R2(K)+2.11E3»Z(M)«R2{K)»EXPF(  1.0065) )/1.E6 

E4  =  2. 1  1E3*Z(M)»EXPF(  1.0065)»(V(I)  +  (2.*R1(J))/1.E3)+V(I)»R1(J) 

F4=8.904E3»Z(M)*Z(M)»EXPF(2.013)+R2(K)»(V( I )-300. ) 

B4=(E4+F4)/1.t3 

C4=4.22»Z(M)*EXPF( 1.0065)»( V{ I )-300. ) 

CT4=(-B4+S0RTF( ( B4*B4 )- { 4. »A4»C4 ) ) )/ (2.*A4) 

IFICT4-1.5)  14, 10, 10 

10  IF(CT4-2.5)  11,11.14 

11  CMAA=MAX1F(CT1,CT2,CT3,CT4) 
CMIN=MIN1F(CT1,CT2,CT3,CT4) 
DELTA=CMAA-CMIN 
IF(DELTA-. 15)  12,12,14 

12  CMAX1=MAX1F(CMAA,CMAX1 ) 
CMIN1=MIN1F(CMIN,CMIN1 ) 

13  CONTINUE 

600  F0RMAT(F6.0,2F 13.0,2F13.5) 

WRITE  OUTPUT  TAPE  4,600,  V  (  l')„ft  1  (  J  ),  R2  (  K  )  ,CMAX  1  ,  CMI  N  1 

14  CONTINUE 

15  CONTINUE 

16  CONTINUE 
END  FILE  4 
PAUSE  1 
END 

END 
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These  maximuv.  m  neon  currents  do  have  value  for  the  analyst 

after  the  design  has  been  chosen,  and  all  that  re     is  the  determina- 
tion of  the  particular  maintaining  voltage  in  any  one  circuit,   This  will 
be  of  immediate  use  provided  that  the  neon  tuba  voltage  is  close  to  one 
of  the  increments  of  voltage  used  in  the  digital  computer  program;  in 
other  cases,  interpolation  is  required. 

The  digital  computer  program  which  is  flow  charted  in  figure  12  and 
represented  in  FORTRAN  statements  in  figui'e  13  has  a  computer  running  time 
of  8  minutes  and  30  seconds.  This  output  data  still  requires  interpret?" 
tion,  evaluation  and  presentation  in  a  useable  manner  before  it  can  be 
considered  as  a  true  design  asset. 
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6.  tion  ata 

Prior  to  an  extensive      ;  retation  of  the  digital  computer 
printout,  it    >  i£  advantageous  to  consider  the  neon  photoconductor 
circuit  .      ,   Let  us  make  the  following  assumptions :  1) 

The      Uning  voltage  of  the  neon  tube  is  fixed  at  60  volts.   2)  The 
known  parameter    has  only  one  value  namely,  2,5.   3)  Temperature  varia- 
tions are  such  that  the  ambient  temperature  will  vary  between  25°  C  and 
55°  C,   With  these  three  assumptions  two  circuits  will  be  considered,  one 
with  compensation  and  one  without  the  compensating  shunt  photoconductor. 
The  printout  for  the  compensated  and  uncompensated  circuits  are  shown  in 
figure  16.  These  two  printouts  are  of  the  same  format  as  the  main  program 
except  that  the  maximum      rnimum  currents  are  now  solely  attributed  to 
temperature  variations. 

A  plot  of  the  design  parameters  R2  versus  Rl,  which  meat  output 
specifics       namely  neon  current  between  1.5  and  2.5  mil  Liatoperes, 
is  constructed  for  both  the  compensated  and  uncompensated  photoconductor 
circuits.  These  plots  are  figures  17  and  18  respectively.  The  uncom- 
pensated photoconductor  circuit  has  a  minimum  value  of  shunt  resistance. 
This  constraint  is  brought  about  by  the  physical  requirement  of  a  suf- 
ficient voltage  drop  across  the  shunt  resistance  to  fire  the  neon  tube. 
The  compensated  photoconductor  circuit  has  a  dark  photoconduetive  re- 
sistance In  series  with  the  shunt  resistance  prior  to  firing  which  is  in 
che  megaohm  region.  Therefore,  only  the  uncompensated  circuii  requires 
a  i      i  value  o£  shunt  resistance.  Assuming  a   firing  voltage  of  135 
volts  which  is  an  observed  average  value,  then  the  following  calcula- 
tions ior  the  minimum  value  of  R2  at  Z ,  ~   2.5  apply: 
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The  points  when  plotted  produce   the  dashed   line   in   figure   18. 
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Ln  each  plot  in  order  to  observe  the 
effect  of  the  ch&Tv  nrent  due  only  to  temperature  variations. 

These  numbers  with      area  represent  the  difference  between  the  maximum 
and  minimum  nson  current  in  sniiliamperes  for  that  particular  set  of  design 
parameters  (Rl  and  R2).   A  comparison  of  the  two  plots  reveals  that  compen- 
sation is  actually  achieved  by  insertion  of  the  shunt  photoconductor.  A 
second  conclusion  that  can  readily  be  made  is  that  for  a  particular  value  of 
Rl  in  the  uncompensated  area  plot  the  variation  in  current  is  constant  and 
independent  of  the  value  of  the  shunt  resistance  R2.  This  constant  value 
is  also  the  limiting  value  which  the  compensated  area  plot  approaches  as 
the  shunt  resistance  becomes  infinite.  Analytically  this  last  result 
could  have  been  foreseen  since  for  large  values  of  shunt  resistance  the 
effect  of  the  compensating  photoconductor  on  the  total  series  resistance 
diminishes  and  eventually  approaches  zero.  This  simple  neon  photoconductor 
circuit  will  be  employed  later  in  this  section  as  an  illustrative  example 
of  best  design  procedure. 

The  digital  computer  printout  of  the  main  program  whether  tt  is  the 
hyperbolic  or  "second  order"  approximation  is  of  the  same  format  (See 
Figure  15  for  printout  format).   For  each  value  of  series  resistance  Rl, 
those  values  of  shunt  resistance  R2  which  mark  the  upper  and  lower  limits 
of  successful  circuit  combinations  in  the  printout  are  plotted  on  graph 
paper  producing  areas.   In  the  case  of  the  hyperbola  only  two  of  these 
plots  are  necessary,  one  for  the  lowest  maintaining  voltage  encountered 
and  one  for  the  highest  maintaining  voltage  sine.        is  always  nega- 

tive.  The  "second  order"  approximation  expression  is  more  involved  and 

\  Y  Z. 
it  cannot  be  mathematically  proven  to  have  a  negative  for  all 

±  V 

design  sets  to  be  considered.   In  this  case  increments  of  maintaining 

voltage  were  used  in  the  digital  computer  program  and  as  many  plots  as 
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volt;  s  a  got      ^n  for 

a  partic  all  Z's  and  all  T's.   In  both 

approximations  the  resulting  areas,  whether  two  at  in  the  case  of  the 
hyperbola  or  many  as  in  the  second  order,  can  now  be  overlapped  to  pro- 
duce a  final  area  which  is  common  to  all  of  the  individual  areas.  This 
final  area  if  it  exists  has  satisfied  all  of  the  restrictions  imposed  on 
temperature,  2,  maintaining  voltage  and  output  neon  current.  The  resultant 
area  for  the  hyperbolic  approximation  is  shown  in  figure  20. 

In  order  to  observe  the  effects  of  more  restrictive  specifications, 
the  printout  data  of  both  the  hyperbola  and  ?t  second  order"  approximations 
was  plotted  using  those  upper  and  lower  limits  of  the  shunt  resistance 
that  yielded  an  output  neon  current  between  1.8  and  2.2  milliamperes. 
After  having  plotted  the  individual  area  plots,  the  overlapping  procedure 
was  employed  to  fLnd  the  common  area.  A   common  area  is  non-existant.   If 
this  specification  of  output  current  had  been  used  vice  1.5  -  2.5  milli- 
amperes,  then  the  absence  of  an  area  indicates  that  there  is  no  design  para- 
meter combination  that  is  satisfactory  for  the  entire  maintaining  voltage 
range . 

The  amount  of  compensation  effected  can  be  seen  by  analysis  of  the 
difference  in  maximum  and  minimum  neon  currants  due  only  to  temperature 
variations.   Refer  to  the  simple  example  mentioned  previously  where  the 
maintaining  voltage  was  fixed  at  60  volts  and  parameter  Z  fixed  at  2.5 
(See  figure  17).  The  mapped  points  in  the  interior  of  the  area  plot  in 
figure  17  represent  the  difference  between  the  maximum  and  minimum  neon 
currents  in  milliamperes  due  solely  to  temperature  variation  for  that 
particular  set  of  design  parameters  (Rl  and  R2).   Inspection  of  the  plot 
reveals  that  the  deviation  due  to  temperature  for  the  same  value  of  R2, 
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n  Rl  ti  -argest  p<        s  lue  m 

area.   For  the  same  value  of  8.1,  the 
deviation  due  to  tempera I     nas  Its  smallest  value  when  R2  takes  on  its 
smallest  possible  value  inside  the  plotted  area.  The  satna  conclusions 
are  true  for  other  similar  plots.  To  see  the  overall  ef£ect  on  the  main 
program  for  the  hyperbola  and  "second  order"  approximations,  the  difference 
in  neon  currents  due  to  temperature  variations  (called  delta  in  the  computer 
program)  was  restricted  to  first  .015  and  then  .010  milliamperes.  The  re- 
sultant area  plots  for  the  hyperbola  without  delta  and  with  delta  equal  to 
.015  are  shown  in  figure  21.    With  the  program  rerun  using  delta  of  .010, 
a  resultant  area  did  not  exist,   A  comparison  of  the  final  area  curves  with 
and  without  delta  reveals  that  the  left  edge  of  the  resultant  areas  which 
employed  the  delta  restriction  are  moved  further  to  the  right,  where  as 
the  right  edge  is  unaltered. 

The  resultant  area  curves,  since  they  are  in  effect  the  common  area 
of  the  individual  areas,,  cannot  be  mapped  for  the  variation  in  neon  currents. 
Any  point  in  the  resultant  area  meets  all  of  the  specifications  that  were 
required  and  is  a  good  design  set.  To  determine  the  best  part  of  the  final 
design  area  we  will  resort  to  inductive  reasoning.   In  inductive  reasoning, 
a  set  of  individual  cases  is  studied  by  the  experimental  method,  and  from 
the  observations  made,  a  general  principle  is  formed.  The  principle  formed 
by  inductive  reasoning  is  reliable  only  whan  all  possible  instances  have 
been  examined.   In  the  individual  areas,  the  mapping  procedure  resulted  in 
the  conclusions  that  the  best  design  set  would  be  the  smallest  value  of  R2 
and  the  largest  value  of  Rl  that  fell  within  the  plotted  area.  Applying  in- 
ductive reasoning  to  this  information,  the  best  set  in  the  resultant  area 
curve  is  also  the  smallest  value  of  R2  and  the  largest  value  of  Rl.  Consider- 
at ion  also  has  to  be  given  to  the  variation  in  the  design  parameters  Rl  and 
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R2,   !■  variant  then  the  apex  of  the  area  would  be  the  best 

design  set.   Rl  and  R2  may  have  tolerances  themselves  of  from  5  to  10  per 
cent  of  their  nominal  value,  hence  some  point  that  could  be  the  center  of 
a  circle  whose  radius  is  the  resistance  tolerance  and  where  the  circle  fell 
entirely  within  the  resultant  area  at  the  lowest  part,  would  be  the  best 
design  set. 

The  resultant  area  curves  of  the  hyperbola  and  second  order  approxima- 
tions are  almost  identical  in  position  and  shape*  Minor  variations  in  shape 
can  be  attributed  to  the  plotting  of  the  individual  computer  incremental  out- 
put data.  As  an  example  of  this  problem  in  plotting  consider  a  design  in- 
crement change  of  2  kohms  in  Rl  where  the  value  of  the  design  parameter  re« 
suited  in  a  successful  circuit  prior  to  the  increment  change,  but  failed 
after  the  increment  change.  Theoretically  then  the  plotted  curve  could  pass 
through  the  first  point  and  it  also  could  pass  through  any  point  to  the  right 
of  this  original  point  up  to  but  not  including  the  value  of  the  original 
point  plus  the  increment.  Here  then  is  another  factor  in  considering  the 
number  of  design  sets  to  be  investigated  versus  computer  time.  This  factor 
should  not  carry  too  much  weight  since  in  the  final  design,  the  designer 
will  pick  a  point  somewhat  removed  from  the  actual  boundary  curve  itself. 
The  resultant  areas  as  obtained  by  the  two  approximations  should  be  the  same 
for  large  values  of  R2,  The  difference  is  noticeable  where  R2  is  280  kohms 
or  less.   In  this  rarge  of  R2  compensation  is  more  effective  and  the  second 
order  area  is  a  more  accurate  representation.  However,  in  the  final  analysis 
the  differences  between  the  two  approximations  prove  to  be  very  minor. 

Returning  to  figure  1?  where  the  voltage  is  fixed  at  60  volts  and  Z  is 
2.5,  it  appears  that  R2 -*—  e*s  at  Rl  =83  kilohms.  The  digital  computer 
printout  data  when  plotted  reveals  that  the  lower  curve  is  in  effect  the 
locus  of  combinations  of  Rl  and  R2  which  result  in  neon  currents  of  1.5 
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mil  Lar  locus  for  curtents  of  2.5 

millamperes.   These  curves  are  actually  only  approximate  locii  of  these 
constant  current  values  since  the  joints  which  were  used  to  construct  the 

curves  were  obtained  by  using         al  values  of  design  parameters. 

2 

The  system  equation  which  was  of  the  form  AI  +  BI  +  C  s  0  was 

solved  explicitly  on  the  digital  computer  for  I?  as  a  function  of  five 
parameteres.   If  the  system  equation  is  rewritten  so  as  to  express  R2  as 
a  function  of  the  remaining  parameteres,  the  result  using  the  hyperbola 
approximation  would  be: 

=   ^  '     '  

Using  a  neon  current  of  2.5  milliamperes,  a  voltage  of  60  volts,  a 
Z  of  2.5  and  a  mean  temperature  of  40°C,  the  value  of  Rl  to  make  the  de- 
nominator of  the  equation  zero  is  35  kohms.  Using  this  equation  similar 
results  could  be  predicted  for  other  maintaining  voltages,  namely  that  in 
the  plot  of  design  parameter  R2  versus  design  parameter  Rl  there  are  values 
of  Rl  where  R2  goes  to  infinity.   Since  the  resultant  area  is  constructed 
from  the  individual  area  plots  we  cannot  expect  to  have  a  closed  resultant 
design  area  plot. 

The  insight  gained  from  the  interpretation  of  the  digital  computer  data 
for  the  circuit  synthesis  problem  revealed  that  the  constructed  boundary  lines 
are  in  essence  locii  of  constant  neon  currents,  namely  the  extremes  of  the 
specified  neon  current  limits  of  1.5  and  2.5  milliamperes.  This  suggests  an 
alternate  technique  for  solving  the  circuit  synthesis  problem.  This  alter- 
nate technique  lends  itself  to  solution  using  an  analogue  computer. 
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7.   Descr  sctronic  Analogue  Computer  Technique 

The  technique  to  be  employed  with  the  analogue  computer  parallels  the 
development  of  the  digital  computer  technique  to  a  certain  extent.   Known 
and  design  parameters  are  as  defined  previously.  The  resultant  circuit 
equation  is  manipulated  to  express  one  design  parameter  explicitly  as  a 
function  of  the  remaining  design  and  known  parameters.   Increments  of  the 
known  parameters  and  end  points  of  the  known  parameters  are  determined  as 
in  the  digital  computer  technique.   The  explicit  circuit  equation  is  then 
converted  into  an  electronic  analogue.   The  electronic  analogue  computer  is 
set  up  so  that  its  equations  have  the  same  mathematical  form  as  the  equation 
for  the  physical  system  under  consideration.   In  the  computer,  voltages  are 
used  to  represent  various  physical  quantities  such  as  velocity,  position,  etc. 
Scaling  factors  are  used  to  convert  the  physical  quantities  into  appropriate 
numerical  values  in  the  computer.  The  mathematical  functions  of  addition, 
multiplication,  etc.,  are  performed  by  electronic  operational  amplifiers 
using  scaled  voltages  as  inputs. 

Large  variations  in  known  parameters  can  be  readily  undertaken  provid- 
ed the  range  in  variation  is  considered  during  the  scaling  process.   Incre- 
mental charges  in  the  known  parameters  are  accomplished  by  merely  varying 
the  coefficient  potentiometers  that  are  associated  with  the  parameter  under 
investigation.   Radical  changes  in  the  known  parameters  may  have  to  be 
handled  by  varying  the  components  of  the  operational  amplifiers,  or  by  the 
insertion  of  one  or  more  additonal  amplifiers. 

A  design  parameter  that  is  implicit  in  the  circuit  equation  is  varied 
by  the  application  of  a  ramp  voltage  input  causing  the  explicit  parameter 
output  voltage  to  be  the  continuous  solution  for  a  fixed  set  of  known  para- 
meters. Utilization  of  an  X-Y  plotter  will  present  a  graphical  representation 
of  a  design  parameter  voltage  as  a  function  of  another  design  parameter  voltage. 
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Each  time  that  e  wi_lj  lt  on 

the  X-Y  plotter.   The  number  o£  known  parameters  that  can  be  varied  on  sny 
one  sheet  o£  the  inet  when  a  loss  of  identity  in  the 

plotted  curves  results.  This  procedure  is  repeated  until  graphical  plots 
are  obtained  for  every  increment  of  each  of  the  known  parameters.  These  plots 
are  antiscaied  on  the  X-Y  plotter  so  that  the  graphical  plots  which  are  ob- 
tained in  this  manner  are  design  parameters  as  used  in  the  original  physical 
system  prior  to  the  analogue  conversion.  Each  plot  has  only  investigated  one 
or  more  variations  in  known  parameters.   If  all  of  the  plots  are  overlayed, 
the  area  which  is  common  to  the  acceptable  area  of  all  the  graphs  is  the  final 
design  area  desired.  Any  point  in  this  final  common  area  meets  all  of  the  re- 
quired specifications. 

Rearrangement  of  the  basic  circuit  equation  employing  the  hyperbola 
approximation  for  the  photoconductors  resulted  in  the  form: 

Q2  .     aSZ'Z  (  £1  t  }c"0ji)-Ji32(  ■  VZ  (  j?j  f  J 

£1  -f-  \cJBrl) 

where  g...    ^  n    #  .*  •* 

The  neon  current  I  is  now  an  end  point  type  of  known  parameter  rather 
than  the  output  quantity  as  in  the  digital  computer  technique.  The  two  neon 
current  values  of  interest  are  1.5  and  2.5  milliamperes.  Physical  quanti- 
ties such  as  ohms,  amperes,  and  maintaining  voltage  in  the  circuit  equation 
are  transformed  into  voltages  for  use  in  the  electronic  analogue  by  appropti- 
ate  scaling.  Known  parameters  of  temperature  and  Z  are  divided  into  in- 
crements as  before.  A  change  in  the  value  of  any  known  parameter  is  accomp- 
lished by  altering  the  individual  potentiometers  that  are  associated  with 
the  parameter  under  investigation.   Since  the  range  of  all  the  known  para- 
meters vary  by  less  than  s.   factor  of  ten,  it  is  possible  to  utilize  one  analo- 
gue computer  circuit  and  account  for  all  of  the  variations  by  suitable  scal- 
ing of  the  physical  parameters.  A  schematic  of  the  analogue  computer 
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amplifier  sequence  is  illustrated  in  figure  22. 

A  step  voltage  ap]  . irgrating  amplifier   F   produces  a 

ramp  voltage  which  corresponds  to  design  parameter  Rl.  This  continuous  vari- 
ation of  design  parameter  Rl,  with  a  fixed  set  of  temperature,  Z,  maintain- 
ing voltage  and  neon  current  will  produce  a  continuous  solution  for  the  ex- 
plicit design  parameter  R2.  Potentiometers  a4,  a5,  a7,  alC,  all  and  al2  are 
held  constant  at  one  value  throughout  the  problem.  These  constant  potentio- 
meter settings  are  determined  by  scaling  physical  quantities  which  are  in- 
variant. The  setting  of  the  variable  potentiometers  at  any  one  instance 
depend  upon  the  known  parameter  that  is  being  investigated.  Table  #2  lists 
the  dependence  of  each  known  parameter  on  its  associated  potentiometers. 

Table  #2 

Variation  in  I  Potentiometers  a3,  a8,  al5,  al6 

Variation  in  T  Potentiometer  a2 

Variation  in  Z  Potentiometer  al 

Variation  in  V  Potentiometers  a9,  al3,  al4 

Figure  23 
Each  of  the  known  parameters  one  at  a  time  is  systematically  varied  by 
changing  their  associated  potentiometers.  The  two  voltages  in  the  electron- 
ic analogue  which  represent  design  parameters  Rl  and  R2  are  automatically 
plotted  by  inserting  these  continually  varying  voltages  as  inputs  to  an 
X-Y  plotter  during  any  one  solution  which  corresponds  to  a  particular  known 
set  of  parameters.  The  output  of  the  X-Y  plotter  is  a  graphical  representa- 
tion of  design  parameter  R2  versus  design  parameter  Rl.  A  typical  output 
curve  for  Z  (photoconductor  resistance  magnitude  multiplying  factor)  corres- 
ponding to  .4  appears  as  in  figure  24.  This  curve  includes  the  effect  of 
temperature,  maintaining  voltage  and  neon  currents.   Solid  lines  represent 
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maintaining  voltages  of  60  volts  and  dashed  lines  represent  maintaining 
voltages  of  100  volts.  The  three  temperatures  investigated  are  25°C,  40°C, 
and  55 °C  which  are  recorded  from  right  to  left  in  groups  of  three  lines. 
For  each  maintaining  voltage,  the  lover  gr<  up  corresponds  to  1.5  milliamperes 
and  the  remaining  group  corresponds  to  2.5  milliamperes.  Similar  curves  for 
other  values  of  Z  appear  in  figures  25,  26  and  27. 

As  was  mentioned  previously  in  section  four,  the  known  parameters, 
temperature  (T),  and  photoconductor  resistance  magnitude  multiplying  fac- 
tor (Z),  were  established  as  incremental  types  of  parameters  since  -  "  '  — 

and  ■        could  not  be  proven  to  be  of  constant  sign  for  the  entire  range 

J 
of  variation.   Examination  of  figures  24  through  27  discloses  that  for  the 

maintaining  voltage  equal  to  60  volts,  the  three  temperature  curves  of  the 

lower  group  (neon  current  equal  to  1.5  milliamperes)  cross  over  each  other. 

This  cross  over  is  indicative  of  a  point  of  inflection  in  -2L~     or  stated 

in  other  terms  the  value  of  sL  changed  sign  at  the  point  of  temperature 

J  T 

cross  over.   By  overlaying  figures  24  and  25  a  similar  effect  would  be  noticed 
for  Z,  namely  that  for  one  value  of  temperature,  maintaining  voltage  and  neon 
current  the  overlayed  curves  corresponding  to  two  different  values  of  Z 

would  cross  over.  This  cross  over  indicates  where  a  point  of  inflection  ex- 

VT 

ists  for  ■?  •""    .   An  overlapping  technique  as  applied  to  the  digital  computer 

areas  is  used  in  finding  the  area  common  to  the  individual  areas.  The  re- 
sultant design  area  is  figure  28.   'since  all  of  the  points  of  inflection  for 
temperature  and  photoconductor  resistance  magnitude  multiplying  factor  lie 
outside  the  resultant  design  area,  these  quantities  are  in  effect  end  point 
type  of  V.nown  parameters  if  the  term  as  used  here  applies  only  to  the  result- 
ant design  area. 

The  design  area  as  obtained  from  the  analogue  computer  in  conjunction 
with  the  X-Y  plotter  meets  the  requirements  of  circuit  performance  with  any 
set  of  known  parameters.  The  Hues  forming  the  boundaries  of  the  area  are 
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exsc;  as  in 

it&l  coi  eter  combination  cannot  be  determin- 

ed  from  this  design  area  as  was  obtained  from  the  digital  computer  data,  since 
there  is  no  additonal  information  available  about  the  resultant  design  area 
upon  which  the  designer  can  base  a  choice.  This  perhaps  is  its  sole  weakness 
as  compared  to  the  digital  computer  technique. 
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8 ,   Com 

The  ,    lem  of  how  t<        a  circuit  to  meet  specifications  while 
using  parameters  which  b.  ges  of  variation  can  be  solved  using 

either  of  the  computer  techniques  described.  The  designer  has  merely  to 
select  design  parameters  from  a  final  area  with  the  assurance  that  the 
circuit  will  meet  the  output  specifications. 

Both  the  digital  and  the  analogue  computer  techniques  provide  adequate 
design  data.   Each  method  has  advantages;  there  is  little  difference  in  the 
amount  of  time  involved  from  start  of  problem  until  solution.   This  time 
naturally  depends  upon  the  designer's  experience  and  familiarity  with  a 
particular  computer.  The  accuracy  is  comparable  in  either  method;  while 
the  digital  computer  is  more  accurate  for  any  individual  calculation,  the 
analogue  has  the  advantage  of  providing  a  continuous  solution  rather  than 
an  incremental  one.  The  analogue  computer,  when  used  with  an  X-Y  plotter 
provides  the  area  curves  directly  without  the  plotting  necessary  to  reduce 
the  digital  data.   However,  the  digital  computer  can  be  easily  programmed 
to  provide  addit  ional  data  useful  in  analysis  of  a  specific,  designed  cir- 
cuit. There  were  no  logic  circuits  employed  in  the  analogue  solution, 
therefore  the  significant  test  for  delta  I  was  not  made  in  this  solution. 

For  the  specific  circuit  used  to  demonstrate  these  techniques,  it 
should  be  noted  that  the  scheme  for  temperature  compensation  did  achieve 
the  desired  compensation  as  reference  to  figures  17  and  18  will  show.  The 
additional  labor  involved  in  employing  the  more  complex  second  order  approxi- 
mation is  not  warranted  since  there  is  no  significant  difference  between  the 
final  design  areas  obtained  with  the  two  approximations.   Figure  29  compares 
the  results  obtained  with  the  two  approximations  using  the  digital  computer 
and  the  hyperbolic  approximation  using  the  analogue  computer. 
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Additional  investigation  should  be  conducted  to  extend  the  study  to 
include  an  analysis  of  the  circuit's  transient  behavior.   A  method  of  in- 
cluding logic  circuitry  in  the  analogue  circuit  would  permit  the  additional 
test  for  delta  I  to  be  conducted.   A  circuit  involving  the  selection  of 
three  design  parameters  presents  an  interesting  problem.   In  this  esse, 
instead  of  acceptable  area  graphs,  the  result  would  be  a  three  dimension- 
al acceptable  volume  plot  of  the  design  parameters. 
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APPENDIX  I 
MAT'IEMATICAL  DERIVATION  LEADING  TO  SECOND  ORDER  EXPRESSION  FOR 
P'lOTOCONDUCTOR  RE  SI  STANCE -NEON  CURRENT  RELATIONSHIP 

To  obtain  a  mathematical  expression  which  approximates  the  photoconductor 
resistance  variation  with  neon  current  shown  in  curve  B,  the  relationship 
is  assumed  to  be  of  the  form: 

A  %f£- 

The  following  points,  selected  from  the  curve  were  used  to  determine  the 
unknown  coefficients  A,  B,  C  and  D. 

I    in  ma.      1.50  1.80  2.20  2,50 

R  in  Kohra  6.71  5.10  3.88  3.37 

A  ( L  7/  )6  $')   r  3  (  ).  s'xt  o~3)  +    (I  (l.71  k  I  6s J  t  (1) 

V4) 


(2) 


r /\    —    s  ?  XV c"    /?  *-  /.  if  X/ 


(1)    -    (2) 


.  if  J  A  .  *  X/6'  ?  3  +   .  5'/  X  /o    C  ^  O     (3)  .  (4) 
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13    I  I   X  J  o 

Substituting  in  (4a)  and  siwplifying^ 

J?.  6"   X  (|,lx    luJ  C-    ~     ^  *"  ^  "  O 

3*  7  s"  X  •    F, 

D  -         .  8  x  /  o  J  e . 
/ 

'W*   X/<^    /?.  f.    /J  x/  </o 

K  p  -     - 

I.  Vt  i    x/c    X     -  I 


To  incorporate  tha  ©fleet  of  tawparatura,  rewrite  in  tfea  following  forms 
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